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1.  Summary 


The  U.S.  Army  Research  Laboratory  (ARL)  Technical  Report  (TR)  ARL-TR-4803  (Jameson, 
May  2009)  describes  the  development  of  a  prototype  fuel  consumption  algorithm  (FCA)  that 
eventually  will  become  a  module  of  the  aircraft  weather  routing  tool.  In  order  to  properly 
compute  fuel  consumption  estimates,  an  accurate  true  airspeed  (TAS)  must  be  input  to  the  FCA. 
In  the  FCA’s  current  configuration,  the  TAS  being  input  is  simply  a  published  value  that  is  valid 
only  under  standard  atmosphere  (SA)  conditions  at  sea  level.  Under  almost  all  potential  flight 
conditions,  the  atmospheric  density  is  less  than  the  SA  sea  level  value  and  the  actual  TAS  is 
correspondingly  greater.  Standard  TAS  values  published  in  aircraft  operations  manuals  can  be 
10-20%  in  error  or  more. 

A  methodology  has  been  developed  to  compute  an  aircraft’s  actual  TAS  from  information  at- 
hand,  in  order  to  use  the  proper  value  in  fuel  consumption  calculations.  The  ambient  flight 
altitude  air  density  is  required  in  order  to  convert  the  indicated  airspeed  (IAS)  to  TAS,  and 
density  is  derived  using  air  temperature  and  pressure.  While  outside  air  temperature  is 
commonly  measured  on  all  manned  aircraft  and  many  unmanned  aircraft  systems  (UAS),  air 
pressure  is  not.  The  methodology  described  herein  works  from  other  infonnation  to  derive  the 
pressure  value;  namely  altimeter  setting  (AS),  surface  elevation  and  temperature,  and  indicated 
altitude. 


2.  Introduction 


Unmanned  Aircraft  Systems  (UAS)  have  become  a  key  component  of  U.S.  military  operations  in 
theaters  abroad,  and  are  finding  their  place  in  border  patrol  and  other  homeland  security  missions 
within  the  Continental  United  States  (CONUS).  UAS  are  likely  to  have  an  increasing  role  in 
reconnaissance,  surveillance,  communication  and  weapons  delivery. 

The  U.S.  Army  Research  Laboratory,  Battlefield  Environment  Division  (ARL,  BED)  has  been 
developing  an  aviation  impact  routing  tool  (AIR)  to  address  the  issue  of  weather  impacts  upon 
various  UAS  missions.  The  AIR  depicts  regions  of  adverse  weather  impacts  (on  a  map 
background  display)  and  then  searches  for  alternate  routes  to  target  areas  and  return  that 
minimize  those  impacts. 

A  related  factor  that  is  to  be  included  in  AIR  involves  UAS  fuel  consumption.  Fuel  consumption 
rates,  remaining  fuel  quantity,  and,  therefore,  allowable  time  on-target  are  vital  issues  that  must 
be  considered  to  safely  and  most  effectively  employ  UAS  in  their  various  mission  roles.  In 
particular,  as  AIR  proposes  alternate  routes  that  minimize  adverse  weather  effects,  additional 
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fuel  requirements  to  fly  those  routes  must  be  taken  into  account.  A  recent  ARL  Technical 
Report  (TR)  (Jameson,  May  2009)  describes  the  development  of  a  prototype  Fuel  Consumption 
Algorithm  (FCA)  that  addresses  these  issues.  The  prototype  FCA  is  being  designed  such  that 
minimal  operator  input  is  required,  while  a  maximum  amount  of  pertinent  flight  information  is 
produced.  UAS  flight  performance  characteristics,  prevailing  winds,  and  mission  route 
definitions  are  incorporated  into  the  FCA.  Using  these  data,  the  algorithm  computes  ground 
speeds  along  the  various  mission  segments  (and  therefore  flying  times),  fuel  consumed,  and  then 
fuel/time  available  to  be  expended  at  a  target  area.  It  is  hoped  that  the  FCA  will  ultimately  prove 
to  be  a  valuable  decision  aid  for  UAS  mission  commanders. 

This  Technical  Note  describes  the  development  of  a  correction  factor  to  be  applied  to  one  of  the 
inputs  to  the  FCA,  the  True  Airspeed  (TAS),  enabling  more  accurate  computations  of  UAS  fuel 
expenditures  and  mission  flight  times. 


3.  Definitions  and  Background 


3,1  Definitions 

1.  TAS:  As  referred  to  in  this  document,  TAS  is  defined  as:  “True  airspeed  is  the  vector 
difference  of  the  velocity  vectors  of  the  aircraft  and  the  air  mass,  both  with  reference  to  the 
earth’s  surface.  When  determining  the  true  airspeed  of  an  aircraft  under  zero  wind 
conditions  and  in  horizontal  flight,  the  true  airspeed  of  the  aircraft  is  equal  to  the  speed  of 
the  aircraft  relative  to  the  earth’s  surface.  When  determining  the  true  airspeed  of  an  aircraft 
under  non-zero  wind  conditions  an  estimation  of  the  wind  speed  vector  is  used.” 

(Wikipedia,3  23  September  2009)  TAS  is  a  function  of  the  thrust  being  developed  by  the 
aircraft’s  engine(s)/propeller(s)  and  (in  an  opposing  sense)  of  the  amount  of  drag  created 
by  its  airframe. 

2.  Station  Pressure  (Pstn):  The  ambient  air  pressure  measured  by  a  “barometer”  at  a  surface 
weather  station.  Pstn  varies  on  a  diurnal  basis  as  the  temperature  varies,  and  is  also 
influenced  by  the  passage  of  high  and  low  pressure  weather  systems.  The  units  of  measure 
of  Pstn  are  the  hecto-Pascal  (hPa)  or  inches  of  mercury  (inHg).  (1.0  inFIg  =  33.86  hPa). 

3.  Air  Density  (AD):  The  relative  “thickness”  of  the  air,  that  is  a  function  of  air  temperature, 
pressure,  and  (to  a  limited  extent)  moisture  content.  The  unit  of  measure  of  “AD”  herein  is 
kilogram  per  cubic  meter  (kg/m  ).  As  air  pressure  increases  with  constant  temperature,  AD 
increases',  as  air  temperature  increases  with  constant  pressure,  air  density  decreases. 


a  Wikipedia,  The  Free  Encyclopedia  is  an  on-line  encyclopedia. 
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4.  Virtual  Temperature  (Tvjrt):  Although  it  might  be  counter-intuitive,  the  moister  the  air  is, 
the  less  dense  it  is.b  Tvjrt  is  used  in  air  density  calculations  as  a  means  of  compensating  for 
the  effects  of  moisture.  Tvjrt  is  a  fictitious  temperature  that  cannot  be  measured,  but  rather 
is  computed  given  the  measured  temperature  and  moisture  content  of  the  air.  Tvirt  is 
always  equal  to  or  slightly  greater  than  the  measured  air  temperature,  being  equal  when  the 
moisture  content  is  zero  and  increasing  as  moisture  content  increases.  As  moisture  content 
increases,  Tvjrt  diverges  above  the  ambient  temperature  and  when  used  to  compute  density, 
the  resulting  density  value  decreases  accordingly.  The  method  of  computing  Tvjrt  will  be 
described  in  a  later  section. 

5.  Moisture  Content  and  Relative  Humidity  (RH):  The  “moisture  content”  of  the  air  refers  to 
the  amount  of  water  vapor  (water  in  the  invisible,  gaseous  state)  that  is  present.  As 
mentioned  earlier,  the  moisture  content  of  the  air  affects  its  density  (with  increasing 
amounts  actually  decreasing  the  density  value).  There  are  several  adequate  measures  of 
moisture  content;  the  measure  that  will  be  employed  herein  is  the  RH.  The  value  of  RH 
varies  from  0%  (no  water  vapor  present)  to  100%  (the  maximum  amount  of  water  vapor  for 
a  given  temperature  is  present). 

6.  Sea  level  Pressure  (Psi):  The  air  pressure  assumed  to  exist  at  sea  level,  given  an  actual  air 
pressure  measured  at  some  other  altitude.  This  is  computed  from  the  Pstn,  the  altitude  of  the 
weather  station  where  it  is  measured,  and  some  estimate  of  the  average  temperature  over 
the  layer  between  the  station  altitude  and  sea  level.  Psi  is  a  reference  pressure  from  which 
altitude  measurements  are  derived. 

7.  Standard  Atmosphere  (SA):  “The  U.S.  SA,  a  series  of  models  that  give  values  for  pressure, 
density,  and  temperature  over  a  range  of  altitudes.”  (Wikipedia,  23  September  2009.) 

There  are  several  SA  parameters  that  are  pertinent  to  this  study. 

At  sea  level  the  SA: 

a.  Air  pressure=  1013.25  hPa; 

b.  Temperature=15  °C,  288  Kelvin; 

c.  Relative  Humidity=0%; 

d.  Density=1.225  kg/nT 

The  SA  Temperature  Lapse  Rate  (the  rate  at  which  temperature  is  assumed  to  decrease 
with  increasing  altitude)  is:  6.5  C/km 


b 


For  reasons  that  are  not  pertinent  to  this  study. 
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8.  IAS:  “Aircraft  display  an  indicated  airspeed  on  an  instrument  called  an  airspeed  indicator. 
Indicated  airspeed  will  differ  from  true  airspeed  whenever  the  aircraft  is  flying  in  air  whose 
density  differs  from  the  density  at  sea  level  and  15  °C.  Indicated  airspeed  is  used  in 
aircraft  operations  as  the  aircraft  stalling  speed  and  structural  limiting  speeds  are  dependent 
on  indicated  airspeed,  irrespective  of  true  airspeed.  However,  proper  navigation  via  dead 
reckoning  (without  constant  ground  reference)  requires  the  use  of  true  airspeed  and  wind 
corrections.”  (Wikipedia,  23  September  2009.) 

9.  Ground  Speed  (GS):  A  quantity  depicting  the  speed  at  which  the  aircraft  is  traveling  in 
relation  to  the  ground  below.  The  GS  is  the  critical  value  that  ultimately  must  be 
determined  in  order  to  accurately  compute  the  time  required  to  fly  each  mission  leg  (and, 
therefore,  the  fuel  expended). 

10.  Altimeter:  An  instrument  on-board  an  aircraft  that  measures  altitude  above  a  reference 
pressure  level  (such  as  the  calculated  sea  level  pressure).  An  altimeter  is  actually  a 
barometer,  but  its  display  is  not  in  units  of  pressure  but  rather  in  units  of  height  above  the 
reference  pressure  level.  An  altimeter  is  mechanically  geared  to  measure  height  assuming 
a  pressure  variation  with  height  according  to  the  SA. 

11.  AS:  “The  value  of  atmospheric  pressure  to  which  the  scale  of  an  aneroid  altimeter  is  set; 
after  United  States  practice,  the  pressure  that  will  indicate  airport  elevation  when  the 
altimeter  is  10  feet  (3  meters)  above  the  runway  (approximately  cockpit  height).” 
(Wikipedia,  29  September  2009.) 

Thus  the  AS  is  a  pressure  value  that  the  pilot  dials  in  to  a  window  on  the  face  of  the 
altimeter;  more  specifically,  it  is  the  reference  pressure  from  which  the  altimeter  will 
measure  and  indicate  its  altitude.  The  AS  is  a  fictitious  pressure,  in  that  it  is  a  pressure 
value  not  measured,  but  rather,  computed  from  the  ambient  air  pressure  at  an  airfield 
(Psfc),c  and  is  a  pressure  that  is  assumed  would  exist  if  a  hole  could  be  dug  from  the  airfield 
down  to  sea  level  and  a  barometer  placed  at  the  bottom  of  the  hole.  This  is  the  same  as  Psi 
except  it  is  given  in  inches  of  mercury  (inHg)  and  computed  from  the  surface  temperature 
rather  than  an  average  temperature.  The  AS  is  also  assumed  to  be  applicable  for  aircraft 
flying  in  the  general  vicinity  of  the  airfield  from  which  it  is  computed  and  it  is  updated 
periodically  as  Psfc  changes 

12.  Flight  Altitude  (FA)  Temperature  (Tfa)  and  Pressure  (Pfa):  Both  temperature  and  pressure 
are  required  for  density  calculations.  Most  aircraft  have  an  outside  air  temperature  (OAT) 
gauge  from  which  the  pilot  can  obtain  the  Tfa.  However,  aircraft  do  not  have  a  means  of 
directly  measuring  Pfa,  and  it  must  be  computed  from  other  known  parameters,  as  will  be 
described  later. 


c  The  method  of  computing  the  AS  from  the  ambient  pressure  will  be  described  in  a  later  section. 
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13.  MATLAB:  The  MATrix  LABoratory  programming  language.  The  prototype  FCA  and  its 
subprograms  (including  the  subprogram  being  described  herein  to  determine  the  correct 
TAS)  were  developed  in  MATLAB. 

3.2  Background 

One  of  the  required  inputs  to  the  FCA  is  the  aircraft’s  TAS.  As  described  in  ARL-TR-4803  and 
in  section  3.1.1.,  the  TAS  vector  is  used  as  one  component  of  a  vector  triangle,  along  with  the 
wind  vector,  to  solve  for  the  third  vector  that  is  the  aircraft’s  GS.  It  is  the  GS  that  must  be 
accurately  determined  in  order  to  correctly  calculate  the  time  of  flight  and  therefore  the  fuel 
consumed  over  the  course  legs  of  a  mission.  To  date,  the  TAS  used  in  the  FCA  has  simply  been 
a  default  value  as  is  typically  published  in  any  aircraft’s  flight  operations  manual. d  However,  as 
will  be  further  explained  below,  under  most  flight  conditions  the  default  TAS  value  is  typically 
too  low.  A  methodology  is  described  herein  to  correct  the  TAS  in  order  to  more  accurately 
compute  fuel  consumption  estimates. 

As  explained  in  the  definitions,  it  is  the  IAS  that  is  measured  by  the  aircraft’s  airspeed  indicator. 
In  almost  all  flight  conditions,  the  IAS  is  not  equivalent  to  the  TAS.  The  airspeed  indicator 
instrument  responds  to  the  differential  between  the  dynamic  pressure  of  the  ram  air  entering  the 
aircraft’s  pitot  tube  (a  tube-shaped  aperture  extending  from  the  leading  edge  of  the  wing)  and  the 
still-air  within  the  aircraft’s  static  pressure  system.  It  is  mechanically  geared  to  read  correctly  (to 
correctly  indicate  the  TAS)  only  when  the  ambient  air  is  at  sea  level  with  the  SA  density  (1.225 
kg/m  ).  At  flight  levels  above  sea  level,  where  in  most  cases  the  air  is  less  dense,6  the  dynamic 
air  pressure  is  reduced  and  the  airspeed  indicator  reads  too  low  (the  IAS  is  lower  than  the  TAS). 
For  some  UAS  and  typical  flight  altitudes  in  the  battlefield,  the  IAS  could  be  as  much  as  10-20% 
in  error.  Also,  the  published,  default  value  of  the  TAS  could  be  in  error  by  the  same  magnitude. 

The  solution  to  the  problem  of  an  inaccurate  TAS  in  the  FCA  is  to  replace  its  default  value  with 
its  actual  value  by  correcting  the  IAS  for  the  ambient  air  density  conditions. 


4.  Approach 


Several  steps  must  be  taken  before  the  correct  TAS  can  be  detennined.  Those  are: 

1 .  Determine  the  atmospheric  pressure  at  the  FA  of  the  UAS. 

2.  Determine  Tvjrt  at  the  FA. 


d  It  is  a  value  that  is  only  applicable  under  SA  conditions  at  sea  level. 

e  With  very  cold  temperatures,  even  at  altitudes  somewhat  above  sea  level,  the  air  might  actually  be  denser  than  under 
standard  sea  level  conditions.  For  most  UAS  flight  conditions,  however,  this  would  not  be  the  case. 
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3.  Utilize  the  atmospheric  pressure  and  Tvjrt  to  detennine  the  AD  at  the  FA. 

4.  Apply  the  AD  as  a  correction  factor  to  the  IAS  to  determine  the  TAS. 


5.  Constants  and  Conversion  Factors 


The  following  is  a  list  of  constants  and  conversion  factors  that  are  required  in  various  equations 
relating  to  TAS  computations.  (These  constants  and  conversion  factors  were  obtained  primarily 
from  (Hess,  1959;  List,  1971). 


mv 

=  18.016 

[  molecular  weight  of  water  vapor  (gm/mole)  ] 

md 

=  28.988 

[  molecular  weight  of  dry  air  (gm/mole)  ] 

eps 

=  0.6215 

[  “epsilon”  -  ratio  of  mv/md  ] 

Le 

=  597.3 

[  Latent  heat  of  vaporization  (cal/gm)  ] 

Rstar 

=  1.98624 

[  Universal  gas  constant  (cal/mole-Kelvins)  ] 

Tsa0 

=  288.16 

[  SA  sea  level  temperature  (Kelvins)  ] 

P  saO  — 

1013.25 

[  SA  sea  level  pressure  (hPa)  ] 

Dsa0 

=  1.225 

[  SA  sea  level  density  (kg/m3)  ] 

F 

J-^so 

=  6.112 

[  Saturation  vapor  pressure  @  0°C  (hPa)  ] 

aO 

=  -0.03414 

[  Constant  required  for  the  pressure  equation  ] 

al 

=  0.190284 

[  Constant  required  for  the  altimetry  equation  ] 

a2 

=  5.255 

[  Inverse  of  “al”,  also  for  the  altimetry  equation  ] 

a3 

=  8.4423  x  10'5 

[  Third  constant  required  for  altimetry  equation  ] 

L 

=  6.5  x  10'3 

[  Temperature  lapse  rate  in  the  SA  (°C/m)  ] 

R 

=  287.3 

[gas  constant  for  dry  air  (ergs/gm-Kelvins)] 

inHg  hPa 

.  =33.86 

[  Conversion  factor  inches  Hg  to  hPa  ] 

ft  meters 

=  0.3048 

[  Conversion  factor  feet  to  meters  ] 

C_K 

=  273.16 

[  Conversion  from  °C  to  Kelvins  ] 

denfac 

=  0.34838 

[  Conversion  factor  in  density  equation  ] 

6.  Pfa  Calculation 


Presently  there  is  no  direct  way  of  measuring  the  ambient  air  pressure  given  the  instruments  on¬ 
board  either  manned  aircraft  or  UAS;  and  yet,  a  value  for  the  pressure  at  flight  altitude  (Pfa)  is 
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required  in  order  to  compute  air  density  and  subsequently  apply  a  correction  factor  to  the  IAS  to 
obtain  TAS.  Two  parameters  that  are  available  from  most  aircraft’s  instrumentation  are  AS  and 
OAT.  Given  these  values,  a  two-step  process  can  be  followed  in  order  to  compute  the  Pfa_f 

1.  According  to  a  National  Weather  Service  (NWS)  Web  site  (NWS,  29  November  2010)  the 
following  equation  can  be  used  to  compute  the  AS  at  an  airfield  location  given  a 
barometer-measured  air  pressure  value  (Psfc,  in  units  of  hPa)  and  the  known  elevation  of 
the  airfield  above  SL  (ZsfC;  in  units  of  meters): 


For  the  purposes  of  this  study,  two  of  the  parameters  in  equation  (1)  are  known  (AS  and 
Zsfc)  and  therefore  it  must  be  solved  for  Psfc. 


First,  substituting  some  of  the  variable  names  for  the  constants: 


AS  —  PSfC  x 


Next,  solving  for  “a3”  in  one  of  the  terms  in  parentheses,  and  substituting: 


8.423  X  10-5  =  a3 


AS  —  PSfC  x 


Rearranging  to  solve  for  Psfc: 
(note  that:  al  =  l/a2) 


ASal  =  Psfcal  x 


1  + 


ASal  =  Psfc  +  (a3  x  Zsfc) 
P?fc=ASal-(a3xZsfc ) 
Psfc  =  [ASal  -  (a3  x  Zsfc)]a2 


(2) 


1  If  in  the  future  Met  data  instruments  that  directly  measure  Pfaare  employed  on  UAS,  the  two-step  process  described  in 
section  6  will  not  be  necessary. 
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2.  Given  Psfc  found  by  equation  (2),  a  second  equation  may  be  used  to  extrapolate  pressure 
from  the  surface  up  to  the  FA  and  find  Pfa.  It  is  (Hess,  1959): 

(Zfa  —  Zsfc )  =  —  x  Tbar  x  loge  (3) 

In  equation  (3),  most  of  the  variable  names  have  been  changed  from  those  in  the  text  book 
reference  in  order  to  confonn  to  variables  listed  in  this  document,  but  the  resulting  value  remains 
the  same.  Unlike  equation  (2)  and  the  AS  equation  from  which  it  is  derived  (which  are  based 
upon  SA  temperature  and  pressure  variations),  equation  (3)  is  based  upon  the  mean  temperature 
within  a  layer  of  the  atmosphere  (Tbar),  or: 

j,  _  (j fa-+Tsfc) 

1  bar  —  2 

The  quantity  (Zfa  -  Zsfc)  is  the  “thickness”  of  the  layer  from  the  surface  up  to  the  flight  altitude. 
Since  the  surface  elevation  is  known,  the  surface  pressure  has  been  determined,  the  value  for  Zfa 
may  be  read  from  the  altimeter,  and  a  value  for  T(-,ar  can  be  computed;  equation  (3)  can  be 
rearranged  to  solve  for  Pfa. 

Substituting  for  -R/g,  let: 

—  =  —  =  -29.29  -»  aO  =  -3.414  x  10“2 

g  a  0 

(Zfa  ~  Zsfc )  —  —  x  Tbar  X  loge  ~~ 


Rearranging  and  solving  for  Pfa: 


l08efe)  = 


aOx(z^a  Zsfc) 
T  bar 


=  eV 

Tbar  ) 

P  sfc 

( 

^a0x(^/a-^s/c) 

—  ^sfc  ^  ^ 

^  Tbar 

(4) 
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7.  Tvirt  Calculation 


The  effects  of  pressure  and  temperature  upon  density  are  significantly  greater  than  that  of 
moisture  content;  and,  aircraft  currently  have  no  means  of  measuring  moisture  content  (including 
RH).  However,  it  is  expected  that  in  the  relatively  near  future  meteorological  (Met) 
instrumentation  packages  may  be  installed  on  some  UAS  that  are  engaged  in  battlefield  missions 
as  well  as  in  homeland  security  operations  that  will  be  measuring  RH.  Consequently,  the  issue 
of  the  effect  of  moisture  content  upon  density  will  be  addressed  in  this  study.  To  do  so,  the 
method  of  computing  Virtual  Temperature  is  described  below: 

A.  Find  the  Saturation  Vapor  Pressure  (Es_Tfa)  at  the  temperature  at  flight  altitude: 

Applicable  equation: 

'°ge  fi?)  =  (=?)  x  (A  -  ?)  (Hess,  1959) 

Let: 

/ mvxLe\  (  1  1  \ 

expo  —  yRstar)  x  (^273  t yaJ 

Use  Tfa  and  solve  to  find  the  value  of  “expo”.  Substitute  in  “expo”. 


loge  (“Y^)  =  exP° 


Solve  for  the  Saturation  Vapor  Pressure 


fEs-Tfa\  _ 

V  Eso  ) 


Es-Tfa  =  Es0  x  eexpo 


B.  Using  Es_Tfa,  find  the  corresponding  Saturation  Mixing  Ratio  (w5_Tfa):  (Hess,  1959) 


ws  Tf 


epsxEs_T fa 

(j3  fa~Es_T  fa) 


C.  Using  w5_Tfa  and  the  ambient  Relative  Humidity  (RH  as  a  decimal  value),  find  the 
ambient  Mixing  Ratio  (w_Tfa):  (Hess,  1959) 

W-Tfa  =  ws-Tfa  x  EH 

D.  Using  w_Tfa  find  the  value  for  Tvirt  (Hess,  1959): 


Tvirt  ~  Tfa  x 


(l+wTfg/eps) 

(i+w_7ya) 
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8.  Density  and  TAS  Calculation 


Find  the  ambient  density  at  the  flight  altitude:  (List,  1971) 

Dfa  =  derifac  x  (-^)  (6) 

Using  the  ambient  density,  the  SA  sea  level  Density,  and  the  IAS,  find  the  TAS:  (Wikipedia,  23 
September  2009) 

TAS  =  IAS  x  (7) 

ADfa 


9.  Example  TAS  Calculation 


9.1  Assumptions 

An  example  TAS  calculation  will  require  that  the  following  assumptions  be  made: 

(1)  The  altitude  that  a  standard  altimeter  would  have  read  during  a  flight  on  the  case-study  day 
can  be  computed  from  a  variation  of  equation  (2). 

Equation  2  incorporates  a  known  ambient  pressure  value  (the  surface  pressure  in  this  case) 
and  a  known  elevation  above  sea  level  (surveyed  for  the  airport  location)  and  is  used  to 
derive  the  AS.  For  the  purposes  of  this  example,  the  equation  will  be  rearranged  such  that 
the  known  AS  and  a  known  flight  altitude  pressure  (see  the  following  paragraph)  will  be 
incorporated  to  derive  the  altitude  reading  given  by  the  instrument. 

(2)  The  UAS  is  flying  at  an  altitude  where  the  actual  ambient  pressure  is  already  known  to  be 
700  hPa. 

This  pressure  is  being  measured  within  the  altimeter  and  converted  to  altitude  by  the 
instrument  (according  to  the  variation  of  equation  (2)  that  is  derived  in  9.2  below). 
However,  this  pressure  is  currently  not  known  to  the  pilot/mission  commander  of  the  UAS. 
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9.2  Altimeter  equation  derivation 


For  the  purposes  of  this  example  only,  the  following  derivation  is  necessary  to  determine  what 
altitude  an  altimeter  would  have  been  reading  had  the  aircraft  been  flying  at  a  pressure  level  of 
700  hPa.  Beginning  with  equation  (2): 

Psfc  =  [ASal  -  (a3  x  Zsfc)]a2 

Replace  Psfc  with  Pfa  and  Zsfc  with  Zfa. 

Pfa  =  [ASal  -  (a3  x  Zfa)]a2 

Rearrange  to  solve  for  Zfa: 

P#  =  [ASal-(a3xZfa)] 

ASa  1  -  Pfl  =  a3  x  Zfa 


-fa 


(^al-P/aa1) 

a3 


(8) 


9.3  Input  Data 


The  following  input  datum  was  valid  at  0500  Mountain  Standard  Time  (MST)  §,  25  January 
2010,  at  the  700  hPa  pressure  level  ( i.e.,  Pfa  =  700  hPa)  above  the  Santa  Teresa,  NM,  NWS 
observation  site  (EPZ),  approximately  32  km  west-northwest  of  the  El  Paso,  TX,  airport  (ELP). 

Temperature  at  700  hPa  (T7oo  =  Tfa):  +1.0  °C 

It  is  assumed  herein  that  the  temperature  at  700  hPa  remained  constant  for  9  hours,  until  1400 
MST. 

The  following  input  data  were  valid  at  1400  MST,  25  January  2010,  derived  from  the  automated 
weather  observing  site  at  ELP.  The  elevation  Above  sea  level  (ASL)  of  ELP  is  3952  ft.  (1205  m 
ASL). 

•  Temperature  at  the  surface  (Tsfc):  +9.4  °C 

•  AS:  1016.5  hPa  (30.02  in  Hg) 

Other  required  input  data. 

•  RH  at  flight  altitude  (two  samples):  0%  and  100% 

•  IAS: . 100  kts 


g  1200  Universal  Time  Coordinate  (UTC). 
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Applying  equation  (8),  the  following  value  was  obtained  for  the  estimated  indicated  altitude: 

Zfa:  9941  ft.  ASL 

9.4  TAS  Computation 

Using  the  input  data  from  section  9.3,  the  following  values  were  obtained  from  the  TAS 
MATLAB  code: 

•  TAS  (100  %  RH): . 117.3  kts 

•  TAS  (0%  RH): . 11 7.1  kts 

The  difference  of  0.2  kts  between  the  TAS  values  is  consistent  with  effect  of  moisture  on  air 
density.  As  explained  in  section  3.2,  the  more  humid  the  air  becomes,  the  less  dense  it  is  and  the 
greater  the  difference  between  the  IAS  (100  kts)  and  the  TAS.  Obviously,  the  effect  of  RH  on 
TAS  computations  is  minimal  and  could  probably  be  neglected  entirely.  However,  since  it  is 
anticipated  that  an  RH  value  may  be  measured  by  instruments  on-board  UAS  in  the  future,  it  is 
being  incorporated  into  these  computations. 

An  on-line  TAS  calculator  that  is  used  by  general  aviation  pilots  was  accessed  to  compare 
against  the  results  obtained  above  (Paragon,  2009).  The  required  input  was: 


•  Indicated  Altitude: . 9941  ft.  ASL 

•  Altimeter  Setting: . 30.02  in  Hg 

•  Outside  Air  Temperature: . +1.0  °C 

•  IAS: . 100  kts 

•  The  resulting  TAS  from  the  calculator  was: . 1 17  kts 


In  this  example  then,  the  output  from  the  TAS  algorithm  compares  closely  with  the  value 
obtained  from  the  on-line  calculator  and  is  assumed  to  be  correct. 
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10.  Conclusion 


A  methodology  is  described  herein  by  which  the  actual  TAS  value  for  a  UAS  can  be  computed. 
Use  of  a  more  precise  TAS  (as  opposed  to  simply  a  default  value)  in  a  UAS  mission  planning 
system  such  as  the  AIR  tool  (with  the  FCA  implemented)  is  key  in  making  accurate  GS  estimates 
that  would  be  anticipated  over  the  course  of  a  flight  (whether  during  the  planning  phase  or  for 
updates  during  mission  execution).  This  in  turn  allows  a  more  realistic  accounting  to  be  made  of 
projected  fuel  expenditures  for  primary  and  alternate  routes,  thus  maximizing  allowable  time  on- 
target  for  the  UAS  and  enhancing  the  overall  safety  of  the  operations. 

TAS  values  derived  using  this  methodology  are  based  upon  prevailing  Met  conditions,  as  well  as 
the  current  aircraft  IAS  and  FA.  The  required  Met  and  flight  parameters  (OAT,  TSfC,  AS,  FA, 
and  IAS)  are  all  readily-obtainable  during  most  (if  not  all)  UAS  operations.  An  example  is  given 
for  which  the  computed  TAS  compared  almost  identically  with  that  found  using  a  commercial 
Web  site.  During  this  research,  a  number  of  other  TAS  computations  were  compared  to  the  Web 
site,  resulting  in  a  mean  difference  of  less  than  one  knot.  In  many  cases  for  which  realistic  Met 
and  flight  parameters  were  input,  a  TAS  value  was  computed  which  differed  from  the  IAS  by  as 
much  as  10-20%.  Assuming  such  variations  could  also  be  expected  from  a  default  TAS  value, 
errors  of  a  similar  magnitude  would  be  projected  into  GS  and  fuel  consumption  computations  as 
well.  Implementation  of  the  TAS  methodology  could  greatly  reduce  or  eliminate  such  errors. 

The  TAS  methodology  was  originally  coded  in  the  MATLAB  programming  language  for  testing. 
It  has  recently  been  converted  to  the  Java  programming  language  and  is  being  integrated  into  the 
FCA.  The  FCA/TAS  combined  algorithm  is  expected  to  be  incorporated  into  the  AIR  planning 
tool  in  the  near  future  to  allow  operational  feasibility  assessments  of  the  alternate  aircraft  routes 
being  proposed. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


AD 

AIR 

ARL 

AS 

ASL 

BED 

CONUS 

ELP 

Es_Tfa 

FA 

FCA 

GS 

hPa 

IAS 

inHg 

kg/m3 

MST 

NWS 

OAT 

Pfa 

Pstn 

Psl 

RH 

SA 


air  density 

aviation  impacts  routing  tool 
U.S.  Army  Research  Laboratory 
altimeter  setting 
above  sea  level 

Battlefield  Environment  Division 
Continental  United  States 
El  Paso,  TX 

saturation  vapor  pressure  at  flight  altitude  temperature 
flight  altitude 

fuel  consumption  algorithm 
ground  speed 
hecto-Pascal 
indicated  airspeed 
inches  of  mercury 
kilogram  per  cubic  meter 
Mountain  Standard  Time 
National  Weather  Service 
outside  air  temperature 
pressure  at  flight  altitude 
station  pressure 
sea  level  pressure 
relative  humidity 
standard  atmosphere 
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TAS 

true  airspeed 

Tfa 

temperature  at  flight  altitude 

TR 

technical  report 

Tsfc 

temperature  at  the  surface 

Tvirt 

virtual  temperature 

UAS 

unmanned  aircraft  systems 

UTC 

Universal  Time  Coordinate 

w_Tfa 

ambient  mixing  ratio  at  the  flight  altitude  temperature 

vv.S'_T  fa 

saturation  mixing  ratio  at  the  flight  altitude  temperature 

Zfa 

height  at  flight  altitude 

Zsfc 

surface  elevation 
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